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1 . 


INTRODUCTION 


'i 


In this report we present various test cases to exemplify 
the use of the three-dimensional code employed to calculate 
supersonic flow over blunt bodies. 

Problem formulation, mathematical framework and the over- 
all program logic are presented in a separate report (ref. 1). 

In section 2 of this report we present the input data and answer 
listings for the first test case which is for a 20° half-angle 
cone at 15° angle of attack and a free-stream Mach number of 7. 
Section 3 of this report involves the test case for a cone- 
ogive-cylinder at 10° angle of attack and a free-stream Mach 
number of 2.86. A complete listing of the code is given in the 
appendix . 


2. FIRST TEST CASE: CONE 

The first test case is for a 20° half-angle cone at 15° 
angle of attack and a free-stream Mach number of 7 as shown 
below. 
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The input data cards required are listed below. Further 
information about the parameters specified on these cards is 
provided in reference 1. Following this list is the input 
data set used for the first test case. 


Card No. 


Format 


Variables 


1 

815 

NSEG : Number of segment points. 

KIND: Flag for kind of segment. 

0 = sphere or circular ogive 

1 = circular cone 

2 = circular cone with flat cut 

2 

8F10.6 

ZSEG : Z - station initiating segment. 

3 

8F10.6 

RSEG : r - coord, initiating segment. 

4 

8F10.6 

DSEG : Distance from centerline to flat 

chord, initiating segment. 

5 

8F10.6 

ASEG: 4>seg ~ Angle between straight 

down and DSEG initiating segment 

5a 

3F10. 6 

ZC: Z at center of longitudinal arc. 
RC : r at center of longitudinal arc. 
RADIUS: Radius of longitudinal arc. 

(Cards 6-11 are read 

in SUB. INPUT) 

6 

3E15.6, 
5X, 15 

XMACH: Mach number 

ALPHA: Angle of attack (degrees) 

GAMMA: Ratio of specific heats 

NREAL : 0 for perfect gas, -1 for real 

gas (pointed cone starting solu- 
tions are generated internally 
for perfect gas option only) . 

7 

3F10.5 

PHIFD : Meridional angle about which 

points are clustered. 

RK : Meridional clustering parameter 

(0 for no clustering) . 

RJ : Radial clustering parameter (0 

for no clustering) . 
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Card No. 


Format 


Variables 


8 


515 


NIT: 


NIPHI: 


NITER: 


No. of points between body and 
shock (max = 20) 

No. of intervals in meridional 
direction (max = 36 ) 

No. of integration steps desired 
(when ZEND is specified set 
NITER to 99999) 


ICONST(4 9 ): 


NCONE : 


Stepsize is computed every 
ICONST (49) iterations (5 is 
nominal ) 

1 for pointed cone solutions, 

2 for all other geometries 


9 


3F10.5 


CONST (9) 
CONST (4) 
CONST (5) 


Courant No. (usually 0.9) 

Radial dissipation constant 
Meridional dissipation constant 


10 


11 


515 


2F10.5 

315 


DISK 1: 
DISK 2: 
TAPE 1: 

TAPE 2 : 


NTDSOS : 


1 reads solution from tape, 2 
writes solution on tape, 3 does 
nothing (logical unit 12) 

1 reads solution from tape, 2 
writes solution on tape, 3 does 
nothing (logical unit 11) 

1 does nothing, 2 stores body 
shape and writes data on tape 
each Z station, 3 writes data 
only (logical unit 9) 

1 does nothing, 2 reads starting 
solution from punched cards, 3 
stores solution on punched cards 
when exiting (logical unit 7). 

If TAPE2 = 1 and DISK 1 and 
DISK 2 = 2 or 3, a pointed cone 
solution will be generated for 
the perfect gas case only 

0 


ZBS : 


ZFLD : 


increment in z for > 
printing shock and 
body variables (ZBS 
> ZEND if not desired) ^ 

increment in z for 
printing field vari- 
ables (ZFLD > ZEND 
if not desired) 


print 
based on 
z station 
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Card No . 


Format 


Variables 


ITPRTB : 


I TP RTF: 


No. of iterations 
for printing shock 
and body variables 
(ITPRTB > NITER if 
not desired) 

No. of iterations 
for printing field 
variables (ITPRTF 
> NITER if not 
desired) J 


print based 
on number 
of itera- 
tions 


NCASE : If > 0, new case follows 


(The following card contains values used in force and moment cal- 
culations or in shifting the origin of the pointed-cone starting 
solution . ) 


12 


5F10.5 

15 


DIAM: 

Length used in calculating 
ence area; usually maximum 
diameter 

refer- 

ALENGT : 

Reference length used in calculat- 
ing moments 

ZREF : 

Moment reference center 


ZCG: 

Center of gravity location 
static margin calculation 

for 

ZSHIFT : 

The value of Z which corresponds 
to the starting cone origin; if 
no shift set = 0 


IFANDM: 



force and moment calculation 
no force and moment calculation 


(If starting solution is to be read from punched cards (TAPE 2 = 
2) , the following three cards are read in main program. If 
solution is read from magnetic storage device, these are not 
required . ) 


12a 

12b 

12c 


5E15.6 

5E15 . 6 

315, 

4E15.6 


XMACH, ALPHA, GAMMA, RK , PHIFD : 
(Defined above) 

R J : (Defined above) 

NIT, NIPHI, NREAL: 

(Defined above) 


P1INF 

RlINF 

VlINF 

GASCON 


free stream pressure,'' 

free stream density, 

free stream velocity, r 

gas constant (1716.0 
for air) J 


real gas 
option only 
(dimen- 
sional ) 
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Card No 


Format 


Variables 


^ n — n ~ 

[If NREAL = -1, gas tables are placed here and will be read in 
SUB.RGAS(523 cards for equilibrium air)] 


(If TAPE 2=2 punch card starting solution is placed here. The 
first card is the Z station of the starting plane and is fol- 
lowed by flow variables at each node.) 


(The following card(s) is used to change the program control 
variables at preselected longitudinal (Z) stations and is read 
in Program MAIN. At least one card is required if no modifica- 
tions are asked for. In this case ZALTER should then be > ZEND) 


13 


F10.5,I2, 
I3,F10. 5, 
12,13 


ZALTER 

NITA 

NIPHEA 

RJA 

RKA 

PHFDA 


STP : 


DISS1 : 
DISS2 : 

NSWCH1 : 


Z station where altering occurs 

New NIT 

New NIPHI 

New RJ 

New RK 

New PHFD 

f O, stepsize determined automati- 
cally 

>0, value of desired constant 
stepsize 

New CONST (4) 

New CONST (5) 

f 0 , new MacCormack 
1 1, old MacCormack 


NSWCH5 : 



no entropy relaxation 
entropy relaxation 


(The following card is used to initialize the force and moment 
calculations, and is read in SUB.COMPUT. This card is needed 
only if IFANDM = 0. *If NCONE = 2 and IFANDM = 0 this card is 
read before the first card 13, otherwise it is read after all 
card 13 ' s . ) 


14 


6F12.8 


FTX 

FTY 

FTZ 

RMTX 

RMTY 

RMTZ 


initial plane forces in the z, r, 
(J) direction 


initial plane moments in the z, r, 
<J) direction 
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Sample 


Input for First Test Case 


0. 0 



0.0 O'. 0 0.0 0.0 0. 0 2_4 

rr j . 

5 1 

I 

8 ;j 11 12 13 15 !6 17 It 11 20 21 22 23 34 25 ■ 77 21 28 30 31 3? 13 34 35 36 37 : 31 40 41 42 43 44 45 it 47 4t 4j 1 51 52 Vj t-4 55 56 57 5i 58 60 SI 62 1 64 65 66 67 61 6) TO 71 77 73 74 71 * 7? 21 * t' 

10001 . 


13 

It 2 J 4 

5 

1 

8 1‘ n :z 13 14 15 IS 17 It 11 70 21 22 23 24 25 « 27 21 78 33 31 32 33 34 M 31 37 36 38 *0 41 4? «3 44 45 46 <7 41 48 SO 51 52 n A 55 56 K 54 58 SC 61 62 63 64 85 66 67 61 IS 19 71 77 73 74 75 76 77 71 7? *:> 

8 0 0 . 



20 18 0.0 0. 90. . 0 .00 . 00 0 01 

|1 2 i 

J 1 

< 

f 1 m 11 0 13 M tt II 17 11 2C 21 22 23 24 25 76 77 1 21 30 31 32 33 W 35 36 37 31 '40 *1 42 43 44 45 *6 ’ 4t 48 SO SI S2 >j 54 55 56 5* 51 60 6t 62 63 64 65 66 61 68 » II 7J 73 M 75 7* 22 7| 74 «: 

1 . 0 



1 . 0 0 . 0 0 .0 0.0 0 12 

rr tt 

5 6 

1 

1S-.I! 13 14 IS IS IT It !1 20 21 • 23 24 25 26 27 21 71 20 31 ' 03 34 3 3 26 J7 31 35 40 41 • 43 44 45 *5 47 41 48 50 51 52 w 54 55 SS S7 5* 51 66 Si 62 W 6* 65 66 6 7 SI SI 70 71 72 TU M ft K T7 71 71 1C. 

1 0 0 0 0 . 


10000. 750 1500 0 11 

Inn 

5 

< 

1 W 11 17 13 M 15 a H '» 70 21 22 23 74 75 26 27 71 28 30 31 32 23 34 3 5 36 37 31 38 40 4M? 43 44 45 *6 47 a 4j 50 51 52 '»> 54 5 5 56 57 56 58 60 61 62 63 64 65 66 6 7 68 61 70 71 7? 73 74 75 75 ’? n 71 tl‘ 


Z‘ 


3 110 10 

|l 2 J 4 

5 6 

1 

8 l v n 12 13 « 15 16 17 « « 21 21 22 73 74 25 26 27 21 21 30 31 32 33 34 35 36 37 31 39 40 41 42 43 44 4S •« <7 « <1 50 51 52 3j i4 55 5* 57 54 SI 60 61 62 63 64 65 66 6 7 66 6 8 7* 7T 72 n 74 75 76 n 71 75 60? 

0. ? 



0 . 0 0 . 0 9 

E 1 • 

5 6 

• 

3 W H 13 M IS 11 17 If 19 20 71 ' 23 24 25 26 27 28 21 30 31 32 23 34 15 36 37 X M 46 41 42 4 3 44 45 *6 4 7 46 «i 50 51 52 5» 14 55 56 57 56 58 U 61 62 63 64 65 66 67 (| 63 TO 71 72 73 >4 75 76 77 71 71 »o' 


o 


13 1500 5 1 8 

fTTTT 

5 S 

1 

8 W n 12 U M ft It 17 n 11 21 21 22 23 24 » 26 77 28 21 30 31 37 M 34 35 36 3 7 36 38 40 41 47 *3 44 45 4$ 47 *• 48 50 51 52 5> 54 55 56 57 56 S3 60 61 62 63 64 65 Sfi 67 61 63 70 n 72 n M 75 76 77 71 71 IC; 

'? 0 . 0 



0 . 0 0 . 0 7 

in - 

5 < 

r 

1 lv 11 13 W 15 IS 17 II 18 20 21 1 21 24 25 26 27 21 21 38 31 12 33 34 35 36 37 M » 4* tl« 4 3 44 45 «€ 47 41 41 50 51 S2 *w 54 55 56 57 56 31 60 SI 67 C 64 *5 «* 67 68 68 70 71 72 73 74 75 76 7» 71 ft JCi 

7. 0 



15.00 1.4 0 6 

E 1 • 

5 1 

i 

8 U 11 1? !3 M » 1* 17 II 20 21 27 23 24 25 26 27 21 21 30 31 23 34 35 36 37 31 33 <0 41 47 4J 44 45 *6 <7 4| 48 50 51 52 3, i4 55 56 57 5* 58 60 61 62 63 64 O 66 67 61 61 TO 71 7? TJ 74 ft 16 77 71 rt «Vi 

0. 0 



0 . 0 r 

E J.JL 

5 8 

r 

8 .-i 11 13 M 15 1* 17 1* If 20 71 22 23 24 25 » 77 21 21 30 31 33 33 34 35 M 31 31 J8 46 41 42 O 44 45 *6 47 « 4| 50 51 52 5j i4 55 56 57 5* 55 60 61 62 63 64 SS 66 67 61 68 70 71 72 >J M ft 76 77 76 79 tZ'- 

0. 0 



0. 0 4 

E J-L 

5 8 

i 

'»■ a m s :i R » n » 7< n n m » n n n n n n n n jj » n 3i » » « <iu u u « a 47 u w m si sn. m ss si n s. » u ii u u h n « i> ii n n n n u m /s k n » ii r.’ 


37 


848839.70884 3 

| 2 J - 

5 6 

7_l 

1 ij n 12 n M « 17 « II 20 71 » JJ 24 25 » 2? a .1 30 31 32 33 )4 35 K I! 31 39 40 41 4! 4J 44 43 .| (7 41 4| SO 51 52 >, 14 S5 SI 17 Si 31 10 till 13 14 IS U It H H H> n 1: TJ W 25 * 71 31 , 

1 . 0 



1 u u u u . 2 

E J < 

5 1 

r 

S I J II 12 13 14 IS '7 "I 11 23 71 22 23 24 ?S 21 71 » a 10 11 33 13 34 IS H 31 31 29 *0 <1 43 13 4. <5,1 12 II <1 SO 51 S3 •>, U 5S 51 47 Si S! 10 II II 1) $4 IS « 17 II II 20 21 IJ 73 74 IS K 7) a a «c: 


2 


l i CARD 1 

il 2 I < 

5 S 

8 

1 w II 1? U 14 15 IS 17 l| 19 20 71 22 23 24 25 ?S 77 2» 79 30 31 32 13 J4 35 IS 37 11 J8 40 41 42 *3 44 45 *6 4 7 46 49 5C SI 52 >, * 55 56 S7 Si 5 1 60 61 62 63 64 65 66 67 68 S 3 70 H 72 73 7< ft * 77 71 71 iai 


OOOOQOOOOJOOOOOOOOOOOOOQOODOQOOOOOOOOOOOOOCDOOOOOOOOOQOOOOOOOOOOOOOOOOOOOOOOOQOO 

1 2 3 4 5 8 7 I 8 10 It 12 13 M 15 *6 t7 II >J JO 21 22 23 24 25 2* 27 71 29 30 3) 32 33 J4 J5 36 37 38 39 40 <l 42 43 44 45 H 47 4< 4J 50 5' 52 53 54 55 56 57 51 53 60 6 ! 62 S3 S4 (5 «6 S’ it S3 70 71 72 73 74 75 76 77 71 71 10 

M M 1 1 II 1 j 1 1 ] I ,i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 I 1 1 1 I 1 > 1 1 1 1 1 ) M 1 I 1 ! 1 I 1 1 1 1 1 1 N 1 1 11 1 1 1 

2 2 2 2 J 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

1 2 3 4 5 8 7 » 3 15 !l '-2 »3 H 15 16 17 II 13 33 21 22 23 24 25 25 27 28 23 30 31 32 33 34 35 36 3 7 31 3 3 40 41 42 4 3 44 4 5 4 6 4 7 l| 43 50 51 52 53 54 55 56 57 51 S3 60 SI 62 S3 64 65 65 67 44 63 70 /I 7? 73 M 75 76 77 78 73 10 

33333333333333333333333333333333333333333333333333333333333333333333333333333332 

44444444444444444444444444444444444444-444444444444444^4444444444444444444444444 

I 2 3 4 5 I 7 I 9 19 1M2 12 !4 15 IS 17 It 19 20 21 22 22 24 25 25 2? 23 29 20 il 32 33 74 25 M y i; 19 <9 <i 12 4 ] 44 4 ; ij 4 : 4g 49 52 5! 52 53 54 55 56 57 «« 55 E0 61 62 S3 64 65 66 57 68 69 70 71 72 71 74 75 76 77 7 | 7 J JO 

35555555555555555555555555555555555555555555555555555555555555555555555555555555 

56BSSS6S6S66S66S6S6S66666S666BS65665S6S66£5GSBB5B5SS6G666366666fe666S666666666SSE 

1 2 7 4 5 i 7 8 9 10 II !2 U 14 15 16 17 II 19 70 ?: 22 22 74 2: 25 2 23 29 30 31 22 33 ii 55 50 57 54 3S 40 4t 42 43 44 <5 46 47 4J 49 5C *2 55 54 55 56 57 SI 59 60 61 62 S3 64 6 5 6 6 67 61 69 70 71 77 73 74 75 76 77 76 >1 Hi 

77777777777777777777777777777777777777777777777777777777777777777777777777777777 

I I f I I I ! 9 1 ! 8 8 8 8 8 8 8 8 U 3 3 8 8 8 8 8 8 3 J 8 3 8 3 S 8 8 3 i 8 3 8 3 3 6 S 8 8 S ! 3 8 3 8 3 8 8 8 3 t 8 3 8 8 8 8 3 8 3 8 i I I I II I ! I 8 

1 2 J 4 S I 7 I 9 10 It 17 13 14 IS IS 17 18 19 25 21 22 25 24 25 2£ 27 28 2 9 3 9 31 32 3 3 54 5 5 3* 3 7 3t 3 2 4 3 41 4 2 4 ] 44 «• 47 4« 49 50 2 T 5? 5 j *4 55 SC 57 5J 5S SO 6! *2 63 64 63 66 17 6t *1 7Q T1 75 71 74 75 76 T7 71 71 Ui 

S3333335333339999999333399333999S939393333S;9S993339 r '13993339999933993393399939 


CtOgE 34/ 


^TANOARD CQRU SOS ' 1 
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A complete output for this test case is now presented. In 
the following, the velocities are made dimensionless with respect 
to the maximum adiabatic velocity, V m = /2y/y-l • Ptcc/ Ptco an< ^ 
the pressure and density for the complete flow field printout are 
normalized with respect to the free stream stagnation conditions. 
Polar coordinates (z,r,<f>) are used with the corresponding velocity 
components (u,v,w) . The circumferential index is k, the radial 
index is j (j = 3 indicates the body surface) . The output con- 
sists of the following sections: 

1. Printout of the input quantities 

2. Printout of the free-stream velocity field and the 
computational mesh. TAU is the radial computational 
variable, XI is the normalized physical radial variable 
running from 0 on the body to 1 at the shock. TXI = 

, TXIT = |4t^T7Tr- ETA is the computational circum- 
° ^ ° V-LAUJ (prr IA \ 3 (DTIL) 

ferential variable, DTIL = 3 (PHI) ~ ' DTILE = 3 (ETA) ' * 

3. Intermediate printout of shock and body variables con- 
trolled by card 11 . 

4. Printout of the flow field at the final z station. 

5. Line-printer plot of the normalized density field at the 
final z station. 

6 . Printout of the shock and body variables at the final z 
station . 

7. The solution reset to the initial z plane using the 
conical property of the converged flow field. Sections 
4-6 are repeated at z-initial. 

8 . Printout of the force and moment calculations. 
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3. 


SECOND TEST CASE: CONE-OGIVE-CYLINDER 



This problem is run as two stacked cases (NCASE > 0). The 
first part of the solution is to obtain a pointed-cone starting 
solution; this solution is then used to start the integration 
down the ogive-cylinder body. The slopes at the cone-ogive 
intersection are matched, thus causing a ZSHIFT^O. This is a 
case which covers most of the options of the code except that 
force and moment calculations are not made. The input cards 
used are now listed and follow the formats displayed in section 2. 
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Sample Input for Second Test Case 





The output listing for this test case is now displayed. 

The variable definitions and normalizations described in section 
2 are, of course, still applicable. The output consists of the 
following sections: 

1. Printout of the input quantities. 

2. Printout of the free-stream velocity field and computa- 
tional mesh. 

3. Intermediate printout of the shock and body variables 
controlled by card 11. 

4. The converged (starting) flow field reset to the initial 
z plane using the conical property of this converged 
solution . 

5. Line-printer plot of the normalized density field at 
the initial z plane for the ogive-cylinder solution. 

6. Shock and body variables at this initial z plane. 

7. Input quantities for the ogive-cylinder case. 

8. Printout of the free-stream velocity field and computa- 
tional mesh for this case. 

9. Flow field, line-printer plot of normalized density, and 
shock and body variables at z-initial (same as items 4, 

5 , and 6 above ) . 

10. Intermediate printout of the shock and body variables for 
the ogive-cylinder solution (controlled by card 11) . 

11. Intermediate printout of the flow field and line-printer 
plot of the normalized density field controlled by card 
11 . 

12. Flow field, line-printer plot of the normalized density 
field, and shock and body variables at the final z 
station of the ogive-cylinder solution. 
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NSEC. KINO 2 1 1 
ZSEG .10000 100.00000 
RSEG .05N0S JM.GM0H7 
DSEG O.OCCOO O.GCOOC 
HSEG 0.00000 0.00000 
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RADIAL MESH DESCRIPTION 
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